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CN . Abstract 

Starting from the semiclassical reduced-action approach to transplanckian scat- 
Z/y • tering by Amati, Veneziano and one of us and from our previous quantum exten- 

, sion of that model, we investigate the S'-matrix expression for inelastic processes 

CN I by extending to this case the tunneling features previously found in the region of 

classical gravitational collapse. The resulting model exhibits some non-unitary S'- 
matrix eigenvalues for impact parameters b < be, & critical value of the order of the 
_ gravitational radius R = 2G^/s, thus showing that some (inelastic) unitarity defect 

Qj • is generally present, and can be studied quantitatively. We find that S'-matrix uni- 

. tarity for 6 < 6c is restored only if the rapidity phase-space parameter y is allowed 

to take values larger than the effective coupling Gs/h itself. Some features of the 
^ resulting unitary model are discussed. 
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1 Introduction 



^ I The ACV eikonal approach to string-gravity at planckian energies [T] has been recently 
■ investigated in the region of classical gravitational collapse. A simplified version of it — 
the reduced-action model of Amati, Veneziano and one of us [2| — has been extensively 
studied at semiclassical level [21 El H], and has been extended by us (CC) to a quantum 
level [H] . The main feature of such a model is the existence of a critical impact parameter 
b = be of the order of the gravitational radius R = 2G^/s, such that, for b < b^, a. classical 
gravitational collapse is expected to occur, while the elastic semiclassical S'-matrix shows 
an exponential suppression driven by the effective coupling a = Gs/h [2j. This suppression 
admits in turn a tunneling interpretation at quantum level [5] , corresponding to a partial 
information recovery, compared to classical information loss. 

The purpose of the present paper is to further study the CC quantum model, in partic- 
ular its extension to inelastic processes in order to see whether the tunneling suppression 
of the elastic channel is possibly compensated by inelastic production thus recovering 
S'-matrix unitarity. 
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The point above is perhaps the key question that the ACV approach is supposed to 
clarify. Indeed, if our S'-matrix model well represents the original string-gravity theory, 
then unitarity is expected irrespective of whether classical collapse may occur for b < be- 
This could be interpreted as full information recovery at quantum level (compared to 
classical information loss) because the suppression of the elastic channel is compensated 
by the inelastic ones. 

Unfortunately, the situation is not a clearcut one, because of the approximations 
involved in the model. On one hand, the reduced-action approach neglects string and 
rescattering corrections which — as argued in [2| — could come in together because of 
the strong-coupling and, eventually, of the short distances involved. Furthermore the 
quantum-extension of [5] is admittedly incomplete because quantum fluctuations involve 
only the transverse-distance dependence of the metric fields, while keeping the classical 
shock-wave space-time dependence as frozen. Finally, our extension of the S'-matrix to 
inelastic processes is based on a weak-coupling procedure which neglects correlations and 
possible bound states, assumptions which could fail in a strong-coupling configuration. 

Indeed, we find eventually that the model shows a unitarity defect for b < be, which 
is dependent on the rapidity phase-space parameter y, in such a way that unitarity is 
recovered in the y ^ oo limit only. This result is interesting because we do have a non- 
trivial unitary model at large y's and all 6's. But it is puzzling also, because it leaves open 
the question of whether, for moderate y, one of the simplifying assumptions above went 
wrong, or whether instead a unitarity defect is a possible feature of quantum gravity in 
the classical collapse region. 

In order to introduce the subject properly, we summarize in sec. |2]both the semiclas- 
sical ACV results for the S'-matrix and the CC quantum extension, by emphasizing its 
tunneling interpretation in the elastic channel. In sec. [3] we derive an improved integral 
representation of the CC tunneling amplitude which is applicable for any values of the 
?/-parameter, and we discuss the role of absorption for the various regimes of the elastic 
amplitude. We start discussing inelastic processes in sec. HI where we provide two classes 
of S'-matrix eigenstates, one corresponding to a weak-field coherent state which exhibits a 
unitarity defect for b < be, and the other with unitary eigenvalues at all 6's, which requires 
a suitably chosen strong-field configuration. The ensuing expectations on the unitarity 
defect around the elastic channel are compared to the direct path-integral evaluation of 
S^S in sec. |5l We find the ^/-dependent unitarity defect mentioned previously, that we 
have quantitatively evaluated at semiclassical level. We also describe the main features of 
the unitary large- 1/ model, by discussing in sec. E] possible hints of further improvements. 

2 The reduced-action approach to gravitational S'- 
matrix 

2.1 The semiclassical ACV results 

The simplified ACV approach [2] to transplanckian scattering is based on two main points. 
Firstly, the gravitational field gfj,u = 'q^y + hfj,^ associated to the high-energy scattering of 
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light particles, reduces to a shock-wave configuration of the form 

h—\ ^ „ = {27rR)a{x)6{x~) , ^++L.-=o = (2vr/?)a(a;)(5(x+) 



la;+=0 



hij = {T:KfQ{x^ X ) I 5ij 



didj 



h{x) 



(la) 
(lb) 



where a, a are longitudinal profile functions, and h{x) = is a scalar field describ- 
ing one emitted-graviton polarization (the other, related to soft graviton radiation, is 
negligible in an axisymmetric configuration). 

Secondly, the high-energy dynamics itself is summarized in the /i-field emission-current 
TC{x) generated by the external sources coupled to the longitudinal fields a and a. Such 
a vertex has been calculated long ago [51 Ej and takes the form 

- V^n = V^aV^a - ViVjaViVja , (2) 

which is the basis for the gravitational effective action pi [9], [10] from which the shock- 
wave solution ([!]) emerges [I]. It is directly coupled to the field h and, indirectly, to the 
external sources s and s in the reduced 2-dimensional action 
A 1„ „_ . (vri?)2 



27iGs 



d 33 as + as VaVa + 



(-(VV)^-2V0- VH) 



(3) 



which is the basic ingredient of the ACV simplified treatment. 

The equations of motion (EOM) induced by ([3]) provide, with proper boundary con- 
ditions, some well-defined effective metric fields a and h. The "on-shell" action A{b,s), 
evaluated on such fields, provides directly the elastic S'-matrix 



5'ei = exp ( -^A{b, s] 



(4) 



Then, it can be shown PQ [2] that the reduced-action above (where R plays the role of 
coupling constant) resums the so-called multi-H diagrams (fig. [I]), contributing a series of 
corrections ~ (-R^/6^)" to the leading eikonal. 



jc = 




+ 



x = b 




+ 



Figure 1: Diagrammatic series of H and multi-H diagrams. 

Furthermore, the S'-matrix ([1]) can be extended to inelastic processes on the basis of 
the same emitted-graviton field h{x). In the eikonal formulation the inelastic S'-matrix is 
approximateljlll described by the coherent state operator 



S = exp ( T-w4,(6, s) ) exp 

( i27ri?ya / d^a^ h{x)n{x] 



n{x) 



\A{x),A\x!)\ = Yb{x-^ 



[a(fc, A;3)e"=-^ + h.c\ = A{x) + A^x) 



(5) 



(6) 



^The coherent state describes uncorrelated emission (apart from momentum conservation [llj). How- 
ever, the eikonal approach based on eq. ^ also predicts [I] correlated particle emission, which is sup- 
pressed by a power of {Gs/h)Y relative to the uncorrelated one, and is not considered here. 
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where the operator Q{x) incorporates both emission and absorption of the /i- fields and Y 
parameterizes the rapidity phase space which is effectively allowed for the production of 
light particles (e.g. gravitons). 

In the following we take the liberty of considering F as a free, possibly large parameter 
which — for a given value of a = Gs/h — measures the longitudinal phase space available. 
This is a viable attitude at large impact parameters b ^ V Gh because the effective 
transverse mass of the light particles is expected to be of order h/b, i.e., much smaller 
than the Planck mass, thus yielding roughly Y ^ 1. On the other hand, we should 
notice that dynamical arguments based on energy conservation [llj and on absorptive 
corrections of eikonal type, consistent with the AGK cutting rules [12], tend to suppress 
the fragmentation region in a 6-dependent way, so as to constrain Y to he O (1) for 
impact parameters in the classical collapse region b = O {R) . However, such arguments 
do not take into account possible dynamical correlations coming from multi-H diagrams, 
as mentioned in footnote [H It is fair to state that a full dynamical understanding of the 
Y parameter is not available yet, and for this reason we shall consider here the full range 
< F < cx). 

In the case of axisymmetric solutions, where a = a(r^), a = a(r^), = 0(r^) it is 
straightforward to see, by using eq. ([2]), that 7i(r^) = (d/dr^)?-^(r^) = —2aa becomes 
proportional to the a, a kinetic term. Therefore, the action can be rewritten in the 
more compact one-dimensional form 

Ur^( a(r2)s(r2) + a{r'')s{r'') - 2pM - —^{1 - p)A , d = '^'^ 



(7) 

where we have introduced the auxiliary field p(r^) 

p = r'{l~ {2nRf^) , h = Air'<j>)= -^{1 - p) (8) 

which incorporates the 0-dependent interaction. The external sources s(r^), s(r^) are 
assumed to be axisymmetric also, and are able to approximately describe the particle- 
particle case by setting 7rs(r^) = 5(r^), 7rs(r^) = 5(r^ — 6^), where the azimuthal averaging 
procedure of ACV is assumed|§ 

The equations of motion, specialized to the case of particles at impact parameter b 
have the form 

^ ■ ^-eir'-b'), (9) 



2np 2ttp 

P=7^2^ir'-b'), p' + - = l {r>b) (10) 
2p2 p 

and show a repulsive "Coulomb" potential in p-space, which acts for r > b and plays an 
important role in the tunneling phenomenon. By replacing the EOM (Q into eq. ([7]), the 
reduced action can be expressed in terms of the p field only, and takes the simple form 

A=-Gs j dr^ (-1(1 - pf - -e{r' - 6^)^ ^ - dr^ L(p, p, r') , (11) 



^The most direct interpretation of this configuration is the scattering of a particle off a ring-shaped 
null matter distribution, which is approximately equivalent to the particle-particle case by azimuthal 



averaging [2]. 
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which is the one we shall consider at quantum level in the following. 

Let us now recall the main features of the classical ACV solutions of eq. ffTUj) . First, we 
set the ACV boundary conditions p{oo) = 1 (matching with the perturbative behaviour), 
and p(0) = 0, where the latter is required by a proper treatment ^ of the = 
boundaryjfl Then, we find the Coulomb-like solution 



= 6^ + -R^(x + sinhxcoshx — Xb — sinh x;, cosh Xb) , (12) 

to be joined with the behaviour p = p{b'^)r'^ for < 6^. With the short-hand notation 

Xb = xip"^)) h = tanhxfc, the continuity of p and p at = 6^ requires the matching 
condition 

p{b^) = tanh Xb = cosh^ , = h(l - tl) , (13) 

which acquires the meaning of criticality equation. 

Indeed, if the impact parameter 6^ exceeds a critical value 6^ = {?)\/?)/2)B? at which 
eq. ffT^ is stationary, there are two real-valued solutions with everywhere regular field, 
one of them matching the iterative solution. On the other hand, for b < be the "regular" 
solutions with p(0) = become complex-valued. 

The action (fTTj) evaluated on the equation of motion becomes 

^ = log(4L2) - log[^ + 1 - - tl) ' = tanhxb) (14) 

and provides directly the 6-dependent eikonal occurring in the elastic S'-matrix, while the 
corresponding /^(r^) ~ 1 — p provides the inelastic coherent state. 

Real-valued solutions for b < be exist but are necessarily irregular, i.e., p(0) > 0. Due 
to the definition of p = r2[l — (27ri?)20], which has the kinematical factor r^, we see that 
such solutions show a singularity of the cj) field of type ~ —p{0)/r^ < 0, so that one can 
check [T] that the metric coefficient hrr must change sign at some value of ~ and is 
singular at r = 0. 

A clearcut interpretation of the (unphysical) real-valued solutions with b < b^ and 
p(0) > is not really available yet. However, we know that in about the same impact 
parameter region classical closed trapped surfaces do exist, as shown in [131 El E]- K is 
therefore tempting to guess that such field configurations of the ACV approach (which 
are singular and should have negligible quantum weight) correspond to classically trapped 
surfaces. In this picture, the complex-valued solutions with p(0) = (which are regular, 
and should have finite quantum weight) would correspond to the tunneling transition from 
the perturbative fields with p{oo) = 1 and positive p to the "un-trapped" configuration 
with p(0) = 0. This suggestion is incorporated in the quantum level, by defining the 
S'-matrix as the path-integral over p- field configurations induced by the action ffTTl) . 

2.2 The quantized CC S'-matrix 

The idea of [S] is to introduce the quantum S'-matrix as a path-integral in p-space of the 
reduced-action exponential. In this "sum over actions" interpretation the semiclassical 



■^A non-vanishing p{0) would correspond to some outgoing flux of V(/) and thus to a (5-function singu- 
larity at the origin of h, which is not required by external sources. 
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limit will automatically agree with the expression in eq. (ITT]) above, which is based on the 
"on-shell" action. Furthermore, calculable quantum corrections will be introduced. 

We thus extend the coherent state definition to the quantum level by introducing 
it in a path-integral formulation where the Lagrangian ( fTTl) occurs, as follows 

p(oo)=l 



where Q{x) acts on the multi-graviton Fock space, but is to be considered as a c-number 
current with respect to the quantum variables p, p. We also assume the ACV boundary 
conditions p(0) = 0, p(oo) = 1 as discussed above. 

In the elastic channel, the fi-dependent exponential in (fT5i) is to be replaced by its 
vacuum expectation value (v.e.v.) 

expj-^y dr . (16) 

Of course, in this quantum extension, no commitment is made to a particular classical 
solution so that the output will presumably contain a weighted superposition of the various 
classical paths satisfying the boundary conditions, that we shall calculate in the following. 
Following the above suggestion, we obtain, in the elastic channel, 

S4b,s) = j ^^^^^^ [Vp{T)] expj-^ j drL,(p,p,r)| . (17) 

p(oo)=l 



where we use the expression (11 II) of the reduced action, with the notations r = r^, 
y = 2Y/ Tc and we introduce the Lagrangian 

Ly=o = L , (18) 

with the boundary conditions p(0) = 0, p(oo) = 1 introduced by ACV and discussed in 
sec. O 

For generic values of y, Ly is complex because of the {1—iy) factor in front of the kinetic 
term, and is thus able to describe absorptive effects due to inelastic production. However, 
in order to deal with a hermitian problem, we start considering the y = limit of 5*61 in 
which Ly is replaced by L, and we shall introduce absorption later on. Although this limit 
for the elastic S'-matrix is somewhat unwarranted — because absorption turns out to be 
very important for unitarity purposes — we shall see in sec. Hlthat the path-integral (1171) 
at y = acquires the meaning of S'-matrix eigenvalue for a class of eigenstates close to 
the vacuum state. Therefore, it is anyway important to discuss it separately. 



Ly{p,P, ^) = 4^ 



;i-iy)(i-p)2-— e(r-62) 
p 



2.3 Elastic S'-matrix as tunneling amplitude 

By then setting y = 0, we shall see that the definition (fT7|l given above is equivalent, by a 
Legendre transform and use of the Trotter formula [15] , to quantize the r-evolution Hamil- 
tonian H{t) to be introduced shortly, and to calculate the evolution operator i/(0, 00), 
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thus reducing the S'-matrix calculation to a known quantum- mechanical problem. In fact, 
by eq. ffTSj) . we can introduce the "conjugate momentum" 

„ dL 1 , ^ , , 

and we obtain 

H{t) ^Up-L = ^ - 1 + ^Qir - b')^ , p = 1 + 2Gn (20) 

from which the classical EOM ffTUl) can be derived. Then, quantizing the evolution ac- 
cording to eq. (fTTl) amounts to assume the canonical commutation relation 



M^m, p.-mc-.---. a^- (21) 



and to quantize the Hamiltonian ( 1201) accordingly: 



h 4a dp^ \ p J h p 



Finally, the path-integral ( fTTj) for the S'-matrix without absorption is related by Trotter's 
formula to a tunneling amplitude involving the time-evolution operator U{0, oo): 

5(6,s) ~r(6,a) = (p = 0|W(0,oo)|n = 0) , Ho\U = 0) = 0, (23) 

where the initial (final) state expresses the boundary condition p(oo) = 1 (p(0) = 0) 
and U{t,oo) is the evolution operator in the Schrodinger picture, calculated with r- 
antiordering. The result ( l23l) expresses the elastic S'-matrix as a quantum mechanical 
amplitude for tunneling from the state |n = 0) at r = oo to the state |p = 0) at r = 0. 

We note that the commutation relation (121 p does not follow from first principles, 
but is simply induced by the path-integral definition (1171) . Note also that here we al- 
low fluctuations in transverse space, but we keep frozen the shock-wave dependence on 
the longitudinal variables x^. This means that our account of quantum fluctuations is 
admittedly incomplete and should be considered only as a first step towards the full quan- 
tum level. This step, defined by (I17p - fp^ . has nevertheless the virtue of reproducing the 
semiclassical result for a — ^ oo. 

A more detailed expression of the tunneling amplitude ( !23l) can be derived by intro- 
ducing the time-dependent wave function 

7/>(p,r) = (p|W(r,oo)|n = 0) (24) 

such that 

T{b, a) = {p = 0|W(0, oo)|n = 0) = -^{0, 0) . (25) 

Since the Hamiltonian ( |20l) is time-dependent, the expression of the wave function at time 
r = is related to the evolution due to the Coulomb Hamiltonian He = Hq + Gs/p by 

iV^(r)) = exp (^^] We(0, oo)|n = 0) (r > b') (26) 



h 

h I V h 



exp [ '^'^^[ exp ( We(0,oo)|H = 0) (r < b') . (27) 
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where, according to eq. (l22l) . we have used "free" evolution for r < 6^. Therefore, the 
tunnehng amphtude is obtained by setting r = in eq. fl7r|) as follows 

r(6, a) = {p = 0|^(0)) = (p = 0| exp exp f^^^ Wc(0, oo) |n = 0) 



■ (28) 

This expression is related, by convolution with the free Gaussian propagator, to the 
function 

^,(p) = (p|We(0,oo)|n = 0), (29) 

which turns out to be a continuum Coulomb wave function with zero energy. In fact, due 
to the infinite evolution from the initial condition 11 = <^=^ p = 1, '?/'c(p) is a solution 
of the stationary Coulomb problem 



1 d^ /I 



4q; dp2 Vp 



+ a - - 1 



Mp) = (30) 



with zero energy eigenvalue (where from now on we express p, r^, 6^ in units of = 4G^s). 
The form of ipcip) is better specified by the Lippman-Schwinger equation 

Mp) = e'*"'' + «Go(0) pvQ^ Mp) , Go{E) = [E - + le]-^ (31) 

and thus contains an incident wave with p = 1, plus a reflected wave for p > and a 
transmitted wave in the p < region. Note the principal value determination of 1/p 
which is important for hermiticity purposes. 

We then conclude that the amplitude (|23l) is, by eq. ( |28l) . the convolution of a gaussian 
propagator with the Coulomb wave function ipdp), which has a tunneling interpretation 
with the Coulomb barrier. In fact, by eq. flHTj) . it contains a transmitted wave in p < 
(where the Coulomb potential is attractive) and incident plus reflected waves in p > 
(where it is repulsive). Calculating ipdp) allows to find an explicit expression for the 
tunneling amplitude (sec. [3]). 

Note that, at 6 = we simply have T(0, a) = V'c(O), so that the tunneling interpreta- 
tion is direct and recalls the well-known problem of penetration of the Coulomb barrier in 
nuclear physics [16]. On the other hand for 6 > 0, the convolution with the free propaga- 
tor changes the problem considerably, and is the source of the critical impact parameter, 
as we shall see below. 



3 Tunneling interpretation and elastic amplitude 

The main purpose of this section is to improve the similar calculation of [5] , by obtaining 
an integral representation of the amplitude which is valid for any values of b and y, even 
the large-?/ region which is important for unitarity purposes (see sec. H]). 

We start calculating the tunnehng amplitude (pSj) without absorption in terms of the 
wave function ( |29|) . We shall then introduce absorption according to the definition (fTSl) . 
by discussing in particular the ^-matrix in the elastic channel. 
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3.1 Basic tunneling wave function 



The explicit solution of flHOj) is given by a particular confluent hypergeometric function of 
z = —Aiap defined as follows 

ipc = NcZ e-^/^$(l + ia, 2, z) , z^" + (2 - z)^' - (1 + ia)$ = 
$ ~ ^-(^+'")(l + 0(l/;z)) , (i2~p->-oo) (32) 

where $ is defined in terms of its asymptotic power behaviour for p — oo and the 
normalization factor to be found below, is chosen so as to have, asymptotically, a 
pure-phase incoming wave for p ~ ^ 1, being an IR parameter used to factorize 
the Coulomb phase. We shall call this prescription as the "Coulomb phase" normalization 
at 6 = oo. 

Here we note that the value c = 2 in $(1 + ia,c, 2;) yields a degenerate case for the 
differential equation in flH^ in which the standard solution with the p — *• — 00 outgoing 
wave, usually called f/(l + ia, 2, z) p2], develops a z = singularity of the form Ajz \ 
B log z. Then, the continuation to p > is determined by requiring the continuity of wave 
function and its flux at p = 0, as is appropriate for the principal part determination of the 
"Coulomb" singularity ( 1311) . The outcome involves therefore an important contribution 
at p > of the regular solution F(l + ia, 2, z), so that we obtain 



ze 



= ze-'/^ (u{l + ia, 2, z) + ^^|^F(1 + ia, 2, z)] (33) 
V r(ia) J 



g(^a-./2) cosh('Ka)z-'" + £Li^e(™+"/2) sinh(7ra)(-2)'" (iz +00) 

rfia) 



We are finally able to determine the normalization factor and the value of ipdO), which 
is finite and non-vanishing, as follows 

r(l-l-ia) r(l + la) coshTra 

a value which is of order e~'^", the same order as the wave transmitted by the barrier. 

3.2 Integral representation of tunneling amplitude at 6 > 



For b > 0, the calculation of T in ( 1281) involves a nontrivial integral, which should be 
investigated with care. A convenient way to perform such calculation uses the momentum 
representation of the Coulomb wave function ipc in which p = t is diagonal. More precisely, 
from eqs. (I19f2ip we introduce the representation {R = 1) 

id ^ i d 

The Fourier transform ipc{t) is defined by 

/+00 
dte'^'^^'Mt) ■ (36) 
-00 

From the stationary Hamiltonian ( l30l) in t-space 

= ah(^f-l + (37) 
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we derive the following differential equation for ipc{t)'- 

dtipcit) 2(ia; — t) ia — 1 ia + 1 



t + 1 



whose general solution is 



^|Jc = N{a) (t-iy^-'it + l) 



\ — ia— 1 



(3^ 



(39) 



In order to have a meaningful integral in eq. (l36l) . we need to shift the singularities 
of (1391) at t = ±1 slightly off the real axis. By shifting both of them upwards, we obtain 
an integral representation for the F-part of ipc'- 



■ r+oo 

^,et(p) = - dt e'2"''* (t - 1 - iOy-\t + 1 - iO)- 

J-oo 



■ia—1 



le" 



TT 



dt 



- 1 



t 



(40) 



where the "retarded" subscript, according to standard Green function notations, indicates 
that the integration contour lies below the singular points of the integrand (as shown in 
fig. [2k), yielding a vanishing result for i;s oc p < 0. The [/-part of ifjc can be obtained with 



-1 



1 



-1 



(a) 



(b) 



Figure 2: Cuts and integration paths for the "retarded" (a) and "causal" (b) solutions of 

eq. (Eg;. 



a "causal" prescription for the pole shift, as shown in fig. [2b: 
^eau(p) = / dt e'^^P' (t - 1 + iO)'"-^(t + 1 - iO)-'"-^ 



— CO 

^2 - 1 llTt 



ze 



-z/2 



oo— iO 

i 



TT 



r(ia) sinh(7ra)f/(l + ia, 2, z) + e-™e(i2)F(l + ia, 2, z) 



(41) 



The Coulomb wave function f[33l) is now easily obtained as a linear combination of the 
retarded and causal solutions: 



ITT 



r(ia) sinh(7ra) 
(4iaL2)i« 



r(ia) sinh(7ra) cosh(7ra 



COSh(7ra)^/'ret - V^cau] 

— cosh(7ra)e'' 



Ct(. 



t2 - 1 V 1 + t 



1 -t 



dt . (42) 
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A convenient representation of ipc with a branch cut at finite t can be obtained by means 
of the relation 



1+tJ \t+l 



and is given bjfl 



^-- TMcosM..) + / )¥^i[tT-i) ^''^ 

It is straightforward at this point to perform the gaussian integration in eq. (1^ 

(vrfe^/ia) V2 

yielding the 6-dependent tunneling amplitude 

^(^ «) = F7T^-r77:r7rrT - / +/ . ^ (46a) 



r(ia) cosh(7ra) \ i-oo-ie A / — 1 
. f .r),H(i^Y e-""-') di (46b) 

r(to) cosh(^a) U_„_„ A jo'l^i + ij " l"""' 

where an integration by part has been performed in the last steplfl 

At = 0, the transition amplitude can be computed by noting that the integral 
I-oo-ie + integrand fl46ap can be closed on the lower half-plane and gives a 

vanishing result. Therefore 

r(0, a) = -J^i^lLy^2 P ('^ - 1)'"'^^° dt = (^^^-^')'" 1 (47) 
r(ia) cosh(7ra) Ji (t + 1)'"+^ r(ia) cosh(7ra) ia 



which correctly reproduces the result in eq. 

On the other hand, at 6 > 0, the integral J^°° is exponentially suppressed with respect 
to _. . This can be shown by bending the paths of the two contributions as shown in 
fig. [3] and by noting that the order of magnitude of the integrand is ~ e~'^" above the cut 
and ~ e^'^" below it. 



3.3 Evaluating absorption at quantum level 

In order to take into account multi-graviton emission, we consider the ^-matrix in eq. ( ITTl) 
with non- vanishing values of the absorption parameter y = 2Y/tt which effectively takes 
into account the longitudinal phase space of gravitons. In the following, we consider y as 
a free parameter [0 < y < oo), independent of a = Gs/h, which can vary from small to 
large values according to the effective transverse mass of the light particles being emitted. 
We note, however, as anticipated in sec. 12.11 that the dynamics (sec. E]) will normally 
introduce correlations, and the latter can depress or emphasize some regions of rapidity 

"^In order to push the integration paths to the point i = 1, a convergence factor (t—iy must be added 
to the integrand whenever the denominator — 1 occurs. 

^This result is exact, and differs eventually by the integration paths from the approximate one in 
eq. (4.19) of [5]. 
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t 



Figure 3: Integration paths of eq. (solid lines). The corresponding deformed paths 
(dashed lines) are such that the lower one dominates T while the contribution stemming 
from the upper one is strongly suppressed. 



phase space, as it happens for the case of energy conservation ItUj, thus providing a- and 
6-dependent constraints on the range of possible y's. 

For y ^ 0, the tunnehng amphtude with absorption T(6, a, y) is again given by eq. ( |25l) . 
but in this case the time-dependent wave function ^IM is determined by a non-unitary 
evolution operator Uy{T, oo), due to the fact that the Hamiltonian operator of the quantum 
system is no longer hermitian, as it should in order to describe absorptive effects due to 
inelastic production. 

In fact, the absorption term in eq. ( |T5i) adds an imaginary part to the kinetic term in 
the Lagrangian fllSp and formally changes the definition of the Hamiltonian and of the 
quantization condition in terms of a new parameter a = a{l — iy): 



H = &{p^-l)+%Q{r-b'), [a/^] = ^, a^a{l-\y). (48) 



A simple way to take into account such changes is to solve the evolution equation for the 
wave-function (t|'?/'(r)) = '?/'(t; r) directly in the momentum representation (1351) in which 
p = t is diagonal. We simply obtain 



i^^(t;r) 



a(r - 1) + ae(r - b^ 



2a dt 



Ht;r). (49) 



For T > b"^, the evolution involves the Coulomb- type Hamiltonian with zero energy 
(due to the boundary condition p(oo) = 1) and we get the solution 

V^(t;r)=(^^j ^^N{a,y), {r > b') , (50) 

where the normalization factor N{a, y) will be fixed later on. On the other hand, for 
T < b"^ we have just free evolution, 

i-^log^{t-r) = &{t'-l) , (51) 



yielding 



i>it; t) = Nia, y) { " _l_e-(i-^)(i-*^)(-^^) , (r < b') (52) 
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and therefore 



t- 1 



1 



t + ll 

By then setting p = and r = 0, we get the desired resuh 

1 



gia(l-i;;)(l-t2)(r-62)gi2a(l-is/)pt 



^(0,0) = iV(a,y) / dt 



t+ 1 



iafe2(i_iy)(t2_i) 



t2 



(53) 



(54) 



which is consistent at y = with the representation f l46ap . and differs from it at y > by 
the replacement 6^ — » 6^ = 6^(1 — iy). 

It remains to determine the proper integration path(s) and the normahzation factor 
in eq. fl54p . In the y = Q limit we require A^(a,0) and the integration path to agree 
with eq. fl46ap . By continuity, the integration path at ?/ > is obtained by rotating the 
original one in counter clock-wise direction, as shown in fig. H^, in such a way to remain 
in the convergence sectors of e'"''^(-^~'^)*^, given by 0/2 < arg(±t) < 0/2 + tx/2 where 
= — arg(l — iy) > 0. 





Re(t) 



Figure 4: (a) Convergence sector and integration paths for the t -representation of the 
tunneling amplitude including absorption, (b) Position of the saddle points fory = 0.5 in 
the complex t-plane. As b approaches infinity, the three saddle points approach the real 
axis at the points —1, and 1. The short red lines indicate the steepest descent directions 
forb= 1/4, 1/2, 1 and2. 

The normalization factor is fixed by the requirement of unitarity at large b, namely 
limfc^oo |'^(&, «,1/)| = 1, and can be determined as follows. Firstly, one notes that the 
integrals along Ci and C2 are dominated by saddle points at ti and t2 respectively, with 
tl — > 1 and t2 — > as 6 — s> 00, as shown in fig. Hh. The saddle point condition is given by 
6^(1 — iy)tk{l — tl) = 1 and one has (cfr. app. A of [5]) 



t+ 1 
2 

k=l 



t2-l 



,iafe2(i_ij^)(t2_l) 



-l)'-'4 



TT 



iatfc(3t| - 1) 



6— +00 



e 



ITT 



ah^il — iy) 



e e 



(55) 
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Secondly, one observes that at large b (and even more at large a) the contribution of the 
saddle point ^2 is suppressed with respect to the contribution from ti, therefore 



T{b,a,y)c^N{a,y) e 



TT 



(6 cx)) . 



(56) 



Finally, from the unitarity requirement, which can be also written as a "Coulomb phase" 
normalization condition 

li- = 1 ^ (57) 

we obtain N{a,y){l — iy)^^" = N{a,0), and we conclude that the elastic ^'-matrix (or, 
the tunneling amplitude including absorption) is given by 

(4iaL2)i°(l-iy)i" /• ei»bHi-iy){t^-i) 
r(ia) cosh(7ra) Jci+C2 t"^ - I \t + 



T{b,a,y) 



dt 



(5^ 



Note that the factor |A^(a, y) \ = e°"'' is needed to cancel the extra large-6 suppression 
and thus enhances the 6 = amplitude e""*^'^"'^-* , which increases to e~™/^ for y ^ oo 









y = 0.1 




y = i 




y = 10 




:-07 1e-06 1e-05 0.0001 0.001 0.01 

Figure 5: Transition amplitude at a = 5 for three values of the absorption parameter 
y. Large-b behaviour in linear scale (a); small-b behaviour in logarithmic scale (b). The 
vertical dashed line in (a) shows the critical value be', the horizontal dashed lines in (b) 
are the boundaries ~ [e~™,e~™/^] of the 6—^0 limits of the amplitude for y ranging from 
zero to infinity. 

In fig. E] we have plotted the dependence on the impact parameter b of the elastic 
S*- matrix, for three values of the inelasticity y. For small j/'s {y = 0.1 say) there are 
some oscillations, due to the interference of the saddle points ti and t2 (on the second 
t-sheet reached across the [—1, 1] cut) for the contour Ci. This shows that the elastic 
unitarity bound is marginally overcome if y is too small. In all other cases (with sizeable 
values of y), we observe that the v.e.v. of S is below 1 thus satisfying the elastic unitarity 
bound, and tends to 1 for large b without oscillations. This is evidence of only one saddle 
point (ti) effectively contributing to the integral for sizeable values of y and b. At larger 
y, fixed b, the vacuum-to-vacuum amplitude is less suppressed than at smaller y, and 
the small-6 suppression of the tunneling amplitude is delayed towards smaller values of 
b < be- Roughly, the turning point is at values of b of order bc{l + y'^)~^^^, thus extending 
to values of b smaller than be the validity of the perturbative behaviour. Nevertheless, in 
the 6 — > limit, the amplitude tends to the (non-perturbative) constant limit e""*^'^"'^-', 
between e~™ {y 0) and e~™/^ (y oo). 
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We thus see the emergence of two absorptive regimes, according to the values of y. 
In the very small-y regime, quantum interference is important, in particular for small 
b — be the saddle points ti and tf^ collide and interfere by confirming the critical role of 
be, but leading to an analytic S-matrix at 6 = be, as explained in |5j. On the other hand, 
for sizeable to large values of y only one saddle point dominates and the perturbative 
and tunneling regimes are hardly distinguishable at 6 ~ be, the perturbative behaviour 
with small absorption being extended to smaller values of b. However, we shall see in the 
following that including inelastic channels will make things even, by restoring the role of 
b = be for unitarity purposes, for any values of y. 



4 Inelastic processes and ^'-matrix eigenstates 

So far we have analyzed the S-matrix in the elastic channel, deriving in eq. (!3H|) an explicit 
expression for the probability amplitude 

r=(0|5|0) (59) 

which represents, in this simplified model (fT5l) of transplanckian scattering, the string- 
string scattering amplitude without graviton emission (a state represented by the graviton 
vacuum |0)). 

We found that starting from the elastic channel (the vacuum state), our quantum 
calculation provides absorption for any value of the impact parameter b, and that for 
b < be (critical value) the tunneling absorption persists even if the graviton-emission 
phase-space parameter y were set to zero. This means that the contribution to the S- 
matrix of quite inelastic states is essential to possibly recover unitarity. 

In this section we investigate the issue of unitarity of our model from various points 
of view. 



4.1 Eigenstates and eigenvalues of the S'-matrix 

A convenient way to determine whether or not the S'-matrix is a unitary operator is to look 
for its eigenvalues. Due to the particularly simple form of our S-matrix as (superposition 
of) coherent state operators in the graviton Fock space, it turns out that the S'-matrix 
eigenstates are functional Fourier transforms of the Fock-space coherent states. In detail, 
we define the generic graviton-coherent-state 

\r]{T)) =e-5(^*''?)exp| / ^ a\x)r]{x'')\ \0) (60) 



where 77 (r) is the distribution function of gravitons in the radial coordinate t = x"^, and 
we have introduced the scalar product notation {rj, () = v{'^)C{'^) cIt. Then, by means 
of a (normalized) functional integration in r-space we introduce the Fourier transform of 
coherent-states 

|{^(^)}) ^ ei(-'-) J mr)] e-'(-'^)|iC) , (61) 
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which are parameterized by the radial function uj{t). It is straightforward to prove 
(app. \M> that such states are eigenstates of the S-matrix 



S\Mr)}) = j [Vp{T)] e-'/^^('''-)^^+'(-'^'')|{^(r)}) = e^"^''^^^''^! |{^(r)}) , (62) 

5p{r) = v/2^(l - p(r)) , (63) 

with eigenvalues e^"^''^'. Furthermore, the a;-states are orthonormal in the continuum 
spectrum and are argued to be complete in the Fock space (app. 1X1) . 

The actual evaluation of the S- matrix eigenvalues involves the path-integral in eq. (lU^ , 
whose action differs from the vacuum one by the cu-dependent contribution (u;,5p). At 
the semiclassical level it is easy to derive the modified equation of motion 

2p- , ^ =-J^dj{T) 64 
p'^ V 

in which io plays the role of external force, depending on the given eigenvalue function 

In the strict ci; = limit we are left with the vacuum state equation characterized by 
the usual matching condition (in the y = limit) 

I = h{i - ti) , [h = m) (65) 

and by pilP') = 4 6^ = pf, = 1/(1 — t^). Real-valued solutions with p(0) = and p{oo) = 1 
exist only for b > be, with b^ = 3-\/3/2. For b < be there are complex solutions, yielding 
a complex- valued semiclassical eigenvalue and a calculable absorption, so that |>S'(a; = 
0;b,a)\ < 1 for b < bc^ This simple observation has the consequence that the S'-matrix 
violates unitarity, to some extent, for values of the impact parameter smaller than the 
critical value be- This means that the y-independent b^ separates the perturbative unitary 
regime [b > be) from a regime where a unitarity defect is possible [b < be), rather than 
separating absorptive and tunneling regimes of the elastic channel, as discussed previously. 
The actual unitarity violation for 6 < 6c is dependent on the relative weight of the small-a; 
states in physical matrix elements and is the subject of the following analysis. 

On the other hand, it is essential to note that, if uj{t) is allowed to take properly 
chosen (large) values, then real- valued solutions of turn out to exist for all b 's, thus 
yielding a real A{u) and a unitary eigenvalue with |5'(co')| = 1. A large class "i?" of such 
solutions is found by setting 



Mr) = fj 



^ (I - p^)Q(r ~ b') + {B - pn)Q{b' - t) 



1 - A 



(66) 



where A G M is arbitrary, pn is the semiclassical solution itself and B = {tb — A) / (1 — A) 
by the continuity requirement on uj and pn a.t t = b"^. By replacing the ansatz (|66l) in the 
equation of motion flM|) we find in the r > b"^ region 

2PR = ^ (r > b') , (67) 
Pr 



^We note that the small-o; solutions with > are singled out by a stability argument so 

that indeed we can have, generally speaking, a unitarity defect and not an overflow. 
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while, for t < b"^, we can take p(r) to be any function with continuous p and p and finite 
p, satisfying p(0) = 0, and matching the Coulomb-hke solution in eq. flU7|) at r = b'^, i.e., 
satisfying p(6^) = tb and p(&^) = p^ = (1 — A)/(l — tl). This is an infinite-parameter set 
of functions, since the Taylor coefficients p^"'\f) (0 < f < 6^) for n > 3 are arbitrary. 

We see that the effect of the parameter A occurring in the external force ujr provided 
by the eigenstate is to renormalize the Coulomb coupling in eq. ( 1671) by the factor 1 — A, 
so that it may become less repulsive for < A < 1 and even attractive for A > 1. 
The main point is, though, that eq. is identically satisfied by the ansatz (lUUj) by 
setting no constraints on p(r) in the < r < 6^ region, so that the external force allows 
automatically real-valued solutions for any value of b. Therefore, for any b, eq. ( I66l) yields 
a family of eigenstates of the S-matrix with unitary eigenvalues depending on an infinite 
set of parameters: two of them (A and tb) characterize the Coulomb problem in eq. fl67j) . 
and an infinity of them (the higher-order Taylor coefficients) span the set of functions 
p(r) for < r < 5^. 

We stress the point that the very existence of such unitary eigenstates is a consequence 
of the quantum structure of the S'-matrix ( |T5l) in which the field p(r) is allowed to fiuctuate 
until it reaches the relevant solution p/j of (!64l) . The only problem of such states {ujr} is 
that their overlap with the vacuum is suppressed by the factor 

|(W|0)p = e-^(™), (68) 

where the exponent is of order a/y. Therefore, such states become important only in the 
y ^ a limit. 

We have thus singled out two families of S'-matrix eigenstates: the small-a; one which 
exhibits a critical value b = be, below which no real-valued semiclassical solutions exist and 
the tunneling phenomenon occurs (with non- unitary eigenvalues), and the large-u; one, 
in which an infinite-parameter family of unitary eigenstates exists, characterized by the 
eigenvalue functions ti;/j(r) in eq. flUUl) . This shows that unitarity is not an exact property 
of our quantum model and indicates that unitarity violations, for any given initial state, 
are determined by the overlap profile of such states on the various eigenstates. 

4.1.1 Sum over eigenstates for the elastic channel 

Using the vacuum wave functional ({cjIIO) = e"^*^"^'"^-* it is easy to construct, by eq. ( j62l) . 
the matrix element 

{0\S\{uj}) = {{uj}\S\0) = e-i("'")e''^(") (69) 
and then, by summing over the complete set |{u;(r)}), the v.e.v. 

{0\S\0) = [ [Vu] {0\S\{u}){{uj}\0) = j [Vuo] e-|(-.-)ei-^(") 

[Vp]e-'^^^P^^^'^^-'^^^''^^''\ (70) 



a result already studied in detail in ref. [5] and in the previous sections. 

We thus remark that the quadratic a;-integration in eq. (ITOll introduces explicitly the 
absorption parameter y in the vacuum equations, via the saddle-point value Us = iSp^r) = 
i^/2ay{l - p(r)). We then recover the equation of motion of the elastic channel 

5(1^ = (71) 

P^ 
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whose solutions are complex for any b value, unlike the uj = limit of eq. ( IMll which 
admits real- valued solutions for 6 > 6c [2] ■ A consequence of this feature is that for any b 
value eq. flTTl) predicts the non-vanishing absorption of sec. 13.31 which, for b < be, tends 
to a finite limit even in the y = limit. Therefore, one has to look in principle at all 
possible inelastic channels in order to check whether the absorption of the elastic one can 
be compensated by the unitarity sum. 



4.1.2 States approximating the unitarity sum 

The simplest approach is to look at the unitarity sum for the 5'-matrix from the point of 
view of the squared matrix elements in eq. (l69l) in order to identify the states that maxi- 
mally contribute to the sum. Since the (quasi)elastic matrix elements are absorbed, and 
the eigenstates with unitary eigenvalues are suppressed by the overlap with the vacuum 
state, the overall unitarity defect is a balance of the two absorptive effects just mentioned. 
A fully quantitative analysis is better done by the method of sec. [5l Here we look at the 
contribution of the unitary eigenstates |{ci^_R}) only and this will provide a lower bound 
to the unitarity sum, as follows 

{0\S^S\0) = [ [Vu] \{{u;}\S\0)\' > [ M |(W|5|0)r = f [Vur] e'^^----) , (72) 



where we have used the fact that [^(tUij)! = 1. Thus the suppression exponent of this 
lower bound is here provided by the vacuum functional {ur,ujr). 

In order to optimize the lower bound above ( 172|) . we look for states that minimize 
{u!r,ujr) in the sample defined by eq. ( l66l) . By imposing stationarity on the infinite set 
of parameters > 3) we easily find that p(r) must be a constant for t < b"^, and 

the latter, by continuity, must be p(&^) = tf,. Therefore, we have the matching condition 

l-A = b\{l-tl) (73) 

which corresponds to a Coulomb problem with "charge" (1 — A). This allows to replace 
A(tf,) in the expression 

1. . a A2 " 



2^'^^''^^^-y(l-A)2 [ 
a 

y 



1 2 

— - tb{l - tf, 



dT + b\i-tt] 

b^ 

(74) 



For b > be, this expression has a vanishing minimum with A = 0, corresponding to 
unitarity fulfillment, with a slope parameter varying from (1 — tf,) ~ 1/(26^) for b ^ b^, 
to tb = tc = l/V^ for b = be (as usual). Instead, for b < b^, the minimum becomes 
non-vanishing, with = it, increasing from tc to U ~ (6^)~^/^ — > oo for b decreasing from 
6c to 0, according to the law 



1 l + Stb-2 



62 2 + 3i 



h'ii + h). (75) 



b 

Correspondingly, the value of A, starting from A = for b = be, increases towards A = 2 
for 6 0, so that the Coulomb potential becomes eventually attractive. The value of (17S1) 
at the minimum becomes 

1.- - ^ 4a (1 - 3t7)2 4a 

-{ujr,ujr) = — ^ — . (76) 

2 y t^{2 + 3tfe)(l + 3tb) y 



18 



and has the property of vanishing in the y ^ oo hmit. 

We tentatively conclude that our quantum S'-matrix is always unitary for b > be and 
may be unitary for b < b^ also, provided y/a ^ oo, the unitarity sum being approximated 
by the urs as given above. This is due to the fact that {ujr, ujr) becomes small in that 
limit, and is consistent with the vanishing of the "unitarity action" Au{y — > oo) that 
we shall derive in the next section. 



5 The unitarity action and its features 

5.1 The unitarity action around the vacuum state 

As an alternative method, it is possible to check unitarity directly by performing the 
sum over S'-matrix eigenstates exactly, at fixed field p(t). Since the integration over u is 
quadratic, the unitarity sum becomes 



{0\S^S\0) = j [Vu] \ {{u}\S\0)\^ = j [Vu] e-|(-'-)e-2^'4('^) 

= J [Vp][V~p] q^I[Up)-Up)] dr-i(5p-5„5p-5,) ^ j [Vp\[Vp\ e'^" , (77) 

where we have performed the tu- integration around the saddle point Us = i{Sp — 6p) = 
iy/2ay{p — p), by introducing the path-integral representation of ^(ti;). It is then straight- 
forward to derive the semiclassical equations 



2p - 2iy{p - p) 

2p + 2iy(p - p) 
which govern the unitarity action 



e(r-6^) 
Q{r-b^) 

p2 



(7J 



Au = - j Hp) - Hp) + iay{p - pf dr . (79) 

From eq. (17HI) we see that, for b > b^ real-valued solutions with p(r) = p(r) exist — 
both equations reducing to the elastic one (ITOl) — for which the on-shell unitarity action 
vanishes, thus implying a unitary S'-matrix, since, at semiclassical level, 

(0|St5|0)semicl = e^^" . (80) 

On the other hand, for b < be, the solutions are necessarily complex and eq. (17S|) can be 
satisfied by setting p = p*, thus yielding the equation 

Q(t - 9) 

2p + %53p = ^ , ^ , (81) 
P^ 

which is equivalent to a coupled set of equations for pi = 9ftp and p2 = ^p. Note that, 
unlike the elastic channel case, the equations flHTl) do not have an analytic structure in p; 
therefore they are to be solved as a coupled set of equations having the form 

= ^^Q{t - b') 

P (82) 

= 3-e(r-62) 
p^ 
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under the boundary conditions 



Pi(0) = P2(0) = pi(oo) - 1 = p2(oo) = . (83) 

We note that the unitarity action (|79l) entering the v.e.v. in eq. (IHOll can be decomposed 
into two pieces: 



iAu = 2 J '^L{p) dr + Aay J dr , {h2 = ^p) . (84) 

The first piece is related to the contribution of the vacuum channel {n = 0) to the unitarity 
sum 

{0\S^S\0) = J2{^\S^n){n\S\0) , (85) 

n 

since, by eqs. (HfTTl) . 

g2/SL(p.) ^ |(o|^|o)|2 = (0|^^|0)(0|^|0) (86) 

where pc is the Coulomb-like solution (fT2|) . The second piece oc h"^ can then be roughly 
interpreted as the contribution to the unitarity sum of the inelastic states, and it will be 
computed in sec. 15.21 

Some simplification in the discussion of ( !82l) is obtained because of the existence of a 
constant of motion of energy type. By multiplying the first equation by p2 and the second 
one by pi and by summing we easily prove the relation (valid for r > b"^) 

^ (('^)' + ^) + = ^ ° ' ^^^^ 

which roughly corresponds to the imaginary part of the single- channel "energy" [pY + 1/ p 
(in the y = limit). 

No additional constant of motion seems to be present, the system appearing to be of 
dissipative type and thus not integrable analytically. We quote a general expression for 
the on-shell unitarity action Au-, derived in app. [Bl 

iX(l/)=2a(^2p2(oo) + 353^^ . (88) 

Here tb = p{h'^) and p2(oo) characterize the given solution, but do not appear to be related 
in closed form, so that no matching condition emerges analytically. Nevertheless, one can 
argue that lAuiy) < with positive y-derivative and that Vmiy^ooAuiy) = 0. Indeed, on 
the basis of the equations of motion one can show (app. IB.ip that 



.dAM 

dy 



4a j pI{t) dr > (89) 



and that, for large y, yp2(T; y) reaches a finite limit R2{t). As a consequence, in eq. flHHl) 
both P2{oo) and t2 = are of order 1/y. It follows that \Au\ = 0{a/y), and thus 
vanishes in the y oo limit. 
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5.2 Numerical results 



We have solved numerically the evolution equations flH^ for {pi, P2), and we have obtained 
the unitarity action (l88l) and the semi classical vacuum-expectation value of S'^S ( ISOl) for 
different values of y. In fig. [6] we show our results for a = 5, y = 0.1, 1, 10, and compare 
them to the elastic quantum S'-matrix squared |(0|S'|0)p = (0|S''''|0)(0|S'|0) which gives the 
vacuum-channel contribution to the unitarity sum (1851) . We shall refer to the solutions for 
|(0|S'|0)p as "exclusive" and to those of (0|S'^S'|0) as "inclusive" over the inelastic states. 

We note that, apart from the unphysical overshoot KOj^lO)! > 1 of the transition 
amplitude at sma\\-y and 6 > 6c0 the inequality I (0|5|0)|2 < |(0|5t5|0)| is always satisfied. 
In the small-y limit, inelastic effects are pretty small, in the sense that | (015*10)1^ ~ 
I (0 1 5'''" 5*10) I . This reflects the fact that pi{T,y) coincide with the vacuum solutions in the 
y ^ limit ffTUj) . Correspondingly, there is a sizeable unitarity violation for b < be — 1.6, 
inelastic effects providing corrections of relative order O [y). 

On the other hand, for large values of y, inelastic effects are very important, and 
the S'-matrix is approximately unitary. In this case, the inclusive solutions are markedly 
different from the exclusive ones. The latter scale as r, y) = p(6^(l — iy), r(l — iy), 0) 
and thus are peaked around r ~ 1/y, with b^ ~ b^/y, as roughly seen in fig. so that 
the tunneling regime is displaced towards smaller values of b. This implies in particular 
that the inelastic weight [cfr. eq. ( !84l) ] y j h\ dr = O (1) thus showing the importance 
of inelastic states yielding a finite (non-vanishing) contribution to the unitarity sum ( |85l) 
in the large- ?/ limit. The inclusive solutions, instead, have hi ~ O {1) around r = 1 and 
h2 ~ 0{l/y) everywhere, yielding a "critical" behaviour around b ~ be, as expected. 

^The small overshoot |(0|5|0)p > 1 at low y = 0.1 for b > 1.5 is due to the oscillations of the quantum 
transition amplitude as seen in fig. [3 compared to the semiclassical evaluation of S'^'S. 
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Figure 6: Comparison of the quantum v.e.v. squared of the S -matrix (dashed lines) with 
the semiclassical v.e.v. of the S'^S operator (solid lines) for a = 5 and various values of 
the absorption parameter y. 
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Since h2 is small for large y values, this implies that the on-shell unitarity action scales 
as a/y, yielding small unitarity violations in this limit (figs. [U|7j) . 

The unitarity action is compared in fig. [7| with the unitarity sum (17211761) provided in 
the previous section. We see that the latter is a good approximation to the unitarity 
action for large y's, thus providing some understanding of the coherent states dominating 
the unitarity sum ( l85l) . with the corresponding inelasticity y. 




Figure 7: Comparison of the unitarity action (solid lines) with the unitarity bound esti- 
mates (7^ ) (dashed lines) for a = 5 and various values of the absorption parameter y. 
On the left we observe that the unitarity violation for b < — 1.6 vanishes for increasing 
values of y. On the right, we see that y times the unitarity action tends to a finite limit, 
which is closely bounded from below by the estimate (76). 



At this point, it becomes important to look at the y oo model, which is unitary. 
Since bl{y) scales as bl{0)/y, unitarity effects are mostly seen in the small-6 region, as 
illustrated in fig. [HI We see that for large y's inelastic effects indeed fill the unitary 
defect. Note that |(0|5|0)| ~ e'^/^ this case (instead of |(0|5|0)| ~ e"™ at y = 0), 
thus showing that inelastic effects compensate a finite unitarity defect around 6 = 0, 
consistently with the previous estimate of y j h"^ dr, providing the order of magnitude of 
such effects. 



6 Discussion 

We have presented here a rather comprehensive study of a quantum extension of the ACV 
gravitational S'-matrix, both for the elastic matrix element (including absorption) and for 
the inelastic ones. We have thus been able to provide an analysis of the unitarity problem 
in the classical collapse region. 

A striking outcome of the paper is that our S'-matrix model satisfies inelastic unitarity 
for all values of b in the large-y limit y ^ a ^ 1. We all know how difficult it is to check 
unitarity, even in well-known theories where no puzzling classical behaviour is present. 
Therefore, this result is a quite non-trivial one and encourages us to further investigate 
the large-y model in detail in order to understand the features of the inelastic production 
which is able to compensate the exponential tunneling suppression in the small-6 region. 

A key role, in recovering unitarity, is played by the quantum structure of our S'-matrix, 
which allows field fiuctuations to build up a class of unitary eigenstates, as explained in 
sec. 14. 1[ Such states, characterized by strong fields and small vacuum overlap at finite 
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?/'s, become actually dominant in the y ^ a limit and turn out to saturate the unitarity 
sum. 

On the other hand, the regime y ^ a = Gs/h appears to be disfavoured for h < he 
on the basis of energy conservation and absorptive corrections pjj, because for b ^ R 
emitted gravitons have a somewhat hard transverse mass ~ h/R, finally restricting y to 
be at most O (log a) and actually O (1) in the classical collapse region]^ This means that 
the unitarity defect that we find for finite ?/'s seems to be the normal feature predicted 
by our model in the physically acceptable range of y^s. An interesting point is that — 
as we noted in sec. [5] — it is a defect and not an overflow. A possible interpretation of 
that would be that, in our quantum model, some information loss does show up in the 
classical collapse region. 

However, we do not really believe the unitarity defect of our model to be a possible 
feature of a consistent quantum gravity theory. We rather think that some of the approx- 
imations of string-gravity theory being used in building up the model were inadequate. 

Perhaps the weakest point of our model is the use of an uncorrelated coherent state 
to represent inelastic production in the S'-matrix for any given field h^r). From the 
original derivation [2], we know that correlations are down by a power of y (actually, a 
power of ay) with respect to uncorrelated emissions. This hierarchy in y could perhaps 
provide a rationale for the need of a large-?/ regime to recover unitarity. Furthermore, 
the existence of correlations could provide a non-linear coherent state, and thus a sort 
of "condensation" field which could change considerably the analysis of saddle-points in 
the strong-field configurations and thus provide a mechanism for recovering unitarity. We 
note that this non-linearity is to some extent predictable from the diagrammatic approach 
of [2] , based on the multi-H diagrams of fig. [H 

We further mention the fact that our quantization procedure keeps frozen the longi- 
tudinal space-time structure of the shock- wave. This also is a weak point, and correcting 
for it — although much more difficult — could provide again further non-linearities in 
the reduced action and in the S'-matrix coherent state. 

A different way of thinking is to believe that — associated to the classically collapsing 
states — there are new quantum states, perhaps bound states, which could contribute to 
the unitarity sum even if the explicit phase-space parameter y were set to zero. We have 
nothing in principle against this point of view, we only find it difficult to implement it in 
a predictive way. 

To sum up, our investigation of the quantum reduced-action model has led, in part, to 
a conclusive answer, by exhibiting a unitary version of the model in the (somewhat formal) 
large-y limit. Future developments include the understanding of the inelastic production 
of the unitary model which is calculable within our approach. Furthermore, in order to 
possibly achieve unitarity at finite values of y, we think we need improvements of the 
model itself, probably in the direction of correlated emission, which looks important at 
finite ?/'s in the classical collapse region. 

^Gravitons (fc) are preferentially emitted in the large-angle region 9k > 9q ~ h/bE {q is the scattered 
particle), so that Y < \og{Eb/fi) if the average graviton number (ug) < 1, or F < \og{Eb/h{ng)) if 
(ug) > 1 (cfr. ref. [H]). By specializing to the collapse region b R, we get the limitation above. 
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A Eigenstates of the ^'-matrix 

In this appendix we determine a set of eigenstates of the quantum S'-matrix, and argue 
that such set is complete in the Fock-space of gravitons. 

The basic ideas are taken from the simpler analogue of a one-dimensional harmonic 
oscillator with destruction and creation operators a and with the usual commutation 
relation [a, a''"] = 1. The bare-bone structure of the 5- matrix ^ is in this case 

S = e'^ , Vl = a + a\ (90) 

where we note that Vl is proportional to the position operator. An eigenvector of 
Vt (and therefore of S) with eigenvalue c<j G M can be formally found by applying to any 
state the operator 5(0. — oj): 

^[5{Sl-u)\ij)]=uj[5{Sl-uj)\ij)] \{uj}) = 6{Q - u})\^) . (91) 

By using the vacuum state {ip) = |0) and the standard integral representation of the Dirac 
delta, we find 

\{uj}) = /" " ^ e-i(^eKia+a^)\Q^ = e-'^^\iO . (92) 



In words, the eigenstates of the position operator can be constructed as Fourier transforms 
of coherent states = e'^*^'^+'^^^|0). In particular, S'ljcj}) = e"^|{co'}). 

It is well known that the set of coherent states \z) : z G C, a\z) = z\z) is (over) 
complete. Actually, also the subset of coherent states involved in eq. flU^ with pure 
imaginary eigenvalues z = i(^ is complete in the Hilbert space H. In fact, the map 
z I— > e^'^^lO) = e'^l^/^l^;), C — if is holomorphic, and thus any coherent state \zq) can 
be represented as a superposition of "pure imaginary" coherent states according to the 
Cauchy integral 

— e^-^'lO) (93) 

27ri Z — Zq 

where z = \(^ runs along the imaginary axis and the sign of e is opposite to the sign of 
'^{zq) in such a way that the integration path can be closed around zq. 

Coming back to the infinite-dimensional Hilbert space of gravitons with the destruction 
and creation operators Aix) and ^^(a:;) in eq. ([H]), we observe that the S'-matrix ffT^ 
involves an azimuthally invariant integration of A{x) + J^{x). It is therefore convenient 
to introduce the canonically normalized operators 

a{r = x^)=r^^ ^ [a{T),a^{r')] = 5{r-r') , (94) 



2v^ VY 

whose eigenstates are coherent states depending on a functional parameter //(r) G C: 

\7]{t)) = e('?'"')-(^*''^)|0) = e-^('''^)e('''"')|0) , a(r)|r/(r')) = r/(r)|r/(r')) (95) 
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with the scalar product notation (rjX) = Jo° v{'^)({'t) dr. We argue, by analogy with 
the one-dimensional case, that the set of coherent states with pure imaginary functional 
parameter //(r) = iCir), ({t) G M, is complete in the Fock space of gravitons. 

With the notations above, the S'-matrix (fTSl) can be written in the compact form 

S= j [Vp{r)] e-/^(^) d.+i(5„a+at) ^ 5^ ^ ^/2^(i _ (gg) 

By using the Backer- Campbell-Hausdorff relations 

for casting operators in normal ordering, we can easily derive the action of the S'-matrix 
on the coherent states: 

5|iC(r)) = / [Vp{T)] e"'/^(^)e-5('^'"'^'')e'('^'""')e^(^'"")e-^('^''^)e'(^'"')|0) 

[Vp{T)] e-'/^(^)e-^K+^P'C+^p)ei(f+'5'"'^')|0) = J [I?p(r)] e-^/^('')|i(C + 5^)) . 

(98) 

In practice, for each path p(t), the coherent state parameter ({t) is shifted by an amount 

In order to look for eigenstates of the S'-matrix, we introduce the functional Fourier 
transform of coherent states 

|Mr)})^iv|[PC(r)]e-^(^'-)|iC(r)), (99) 

where is a normalization factor which can be determined by computing 

,-i(<^',C')+i(^,C)-|(C'-C,C'-C) 



(WmMr)}) = N'*N J [VC{r)][Var)] e 

= Ar'*Are-^('"''"') J [VCir]] e'^^''^"'^') = \N\h-'2(^'^^5{{iu - u'}) (100) 

thus requiring = e^^'^'^) for (r)}) to be a complete and orthonormal set. 

This set diagonalizes the S'-matrix operator. In fact, by using eqs. fl98|99p we find 



S\Mt)}) = N I [Vpir)][Var)] e-'/^(^)e-'(^'-)|i(C + 5p)) 

[Vp{T)] e-'/^(^)+H5p'-) |{cu(r)}) (101) 



where we have decoupled the two integrations by shifting C ^ C = C + ^p- The eigenvalue 
of the S'-matrix relative to the eigenstate |{ct;(r)}) is expressed by a path-integral in p 

_ / p-i/L(p)+i{5p,w) 



eigenvj^) = e'^^^^ = J [Vp{r)] e 



which can be estimated in the semiclassical approximation by finding the path Pt^('r) 
around which the "action" A[uj] is stationary, as explained in sec. (14. ip . 



^This motivates the notation Sp in the definition 
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B The unitarity action 



In this section we compute the unitarity action (ff9l) corresponding to the stationary/ classical 
trajectory determined, for b < be, by the equation of motion ( |82l) and boundary condi- 
tions ( |83l) . In terms of the real components (pi,p2) defined by p = pi + ip2 = p*, the 
unitarity action reads 

POO 

Au = -2ia / (2pip2 - 2p2 + 2ypl - K) dr (102a) 
Jo 

K(pi, P2; r) ^ 0(r - b')Q- = 0(r - 6')^^—^ • (102b) 

P Pl+ P2 

In the interval < r < 6^, the potential T4 vanishes. Therefore, the equation of 
motions Pi = p2 = determine a free evolution for the p field, whose solution is Pfc(r) = 
tfcT, [k = 1,2), where the = Pk{^) are free parameters (eventually constrained by the 
boundary conditions at r = 00), having taken into account the initial condition Pfc(O) = 0. 
The corresponding contribution to the action amounts to 

Au\r<V2 = -Aiab%t2 -t2 + ytj] . (103) 

In the interval t > b"^ the evolution is nontrivial, and we need some relations among 
the pfc's and their r-derivatives. Since the "unitarity lagrangian" in eq. (11021) is time- 
independent for r > 6^, the corresponding hamiltonian 

= 2i[2pip2 + 2ypl + K] = (104) 

is a constant of motion, and evaluates to zero because of the boundary condition f){oo) = 1 
that implies pi(oo) = 1, P2(oo) = 0, Vu{oo) = 0. Another useful relation is obtained by 
multiplying the first equation of (IH^ by p2 and the second one by pi, yielding 

2piP2 + 2pip2 + %pip2 = 3?4^P + ^\^P = ^- = Vu. (105) 

p^ p^ p 

In turn, by using the identities (P1P2) = P1P2 + P1P2 + 2pip2, (P2) = 2p2P2 + 2p| and the 
integral of motion fll04p . we obtain 

2(pip2 + yplY + K = . (106) 

The action for t > b"^ can now be computed: 

Au\r>b^ ^ -2i« / (-2p2 - 2K) dr 
IT061 



b2 



POO 

4ia / [p2 - 2(pip2 + ypl)"] dr 
= -4ia[p2(&') - P2(oo) + 2(pip2 + yplYioo) - 2(pip2 + yp^) '(&')] • (107) 

The values of Pfc(&^) and of its derivatives are matched with those of the free solution 

for r < At r ^ 00 we have pi = 0(r), p2 = C(l), p2 2p2/p? = 0{t-^), 

p2 = {t'^), hence (pip2 + ypD' P2(oo). 
By summing the results fll03lll07p we obtain 



Au = -4ia[p2(oo) - 36^t2(ti + 1/^2)] = -4ia 



. X 3 t2 

P'^^^ - 2^, 



(10^ 
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where in the last equahty we exploited the relation 



Vu{hl) 



h\tl + tl) ■ 



(109) 



obtained from the r — > 6^ limit of the integral of motion (11041) . 
B.l y — ^ oo limit 

The boundary problem defined in eqs. ( I82f83|) admits a well defined limit for y ^ oo. In 
fact, by setting R2{j) = yp2{j), we obtain 



The above system has a finite solution for the pair of functions (pi,i?2) in the y oo 
limit. We deduce that, at large y, the real part pi of p tends to a finite limit, whereas the 
imaginary part p2 of p uniformly scales as ^i?2^~°°'. Therefore, the quantities p2{oo), t2 
and Au linearly vanishes with 1/y. The fact that limy^oo-^« = suggests the unitarity 
of the model at |/ = oo. 
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